Comparison of Some Exact and Perturbative Results for a Supersymmetric 

Gauge Theory 



SU(AQ 



(N 

o 

IX, 



43 ■ 



>: 

On 1 
(N 
t— i ; 

o: 



X 



Thomas A. Ryttov° and Robert Shrock b 
(a) Department of Physics, Harvard University, Cambridge, MA 02138 and 
(b) C.N. Yang Institute for Theoretical Physics, Stony Brook University, Stony Brook, NY 11794 

We consider vectorial, asymptotically free M = 1 supersymmetric SXJ(N C ) gauge theories with 
Nf copies of massless chiral superfields in various representations and study how perturbative pre- 
dictions for the lower boundary of the infrared conformal phase, as a function of Nf , compare with 
exact results. We make use of two-loop and three-loop calculations of the beta function and anoma- 
lous dimension of the quadratic chiral superfield operator product for this purpose. The specific 
chiral superfield contents that we consider are Nf copies of (i) F + F, (ii) Adj, (iii) S2 + S2, and 
(iv) A2 + A2, where F, Adj, S2, and A2 denote, respectively, the fundamental, adjoint, and sym- 
metric and antisymmetric rank-2 tensor representations. We find that perturbative results slightly 
overestimate the value of Nf^ cr relative to the respective exact results for these representations, i.e., 
slightly underestimate the interval in Nf for which the theory has infrared conformal behavior. Our 
results provide a measure of how closely perturbative calculations reproduce exact results for these 
theories. 



I. INTRODUCTION 



A longstanding question in gauge theories concerns 
how well the beta function, calculated to same order in a 
perturbation expansion in the gauge coupling, describes 
the properties of an asymptotically free gauge theory, in 
particular, its evolution from high Euclidean momentum 
scales fi in the ultraviolet (UV) to low ^ in the infrared 
(IR). Part of the difficulty in answering this question 
stems from the fact that only the first two terms (i.e., 
the one-loop and two-loop terms) in the beta function 
are scheme- independent. Recall that the beta function 
effectively resums the naive perturbation expansion in 
defining a running gauge coupling <?(//). There is special 
interest in the case where the two-loop beta function has 
a zero away from the origin. If this zero occurs at a very 
small value of a(fi) = g{n) 2 /(Air), then one expects that 
the theory does not confine or produce bilinear fermion 
condensates and associated spontaneous chiral symme- 
try breaking (SxSB). In contrast, if the theory has suffi- 
ciently few fermions, then this IR zero of the beta func- 
tion occurs at sufficiently large a that one expects the 
theory to confine and spontaneously break chiral sym- 
metry. For ordinary non-supersymmetric SU(iV c ) gauge 
theories, based on the calculations of the one-loop [1| and 
two-loop [H terms in the beta function, there have been 
a number of studies of this UV to IR evolution (an early 
work is [|[ ) . A particularly interesting possibility is that 
the IR zero of the beta function could occur at a value 
only slightly larger than the critical value for SxSB, so 
that the theory would remain quasi-conformal for a large 
interval in ln/x, with a large but slowly running ("walk- 
ing" ) coupling and an associated large anomalous dimen- 
sion 7 m of the bilinear fermion operator [H @ ■ There 
has been intensive recent work to study quasi-conformal 
behavior using lattice methods, and much progress has 
been made [6]. These studies have considered not only 
fermions in the fundamental representation of the gauge 



group (usually SU(2) or SU(3)), but also fermions in 
higher-dimensional representations. This has reflected 
interest in quasi-conformal behavior in gauge theories 
with fermions in higher-dimensional representations; for 
a review, see Q- 

In this paper we consider vectorial, asymptotically free 
Af = 1 supersymmetric SU(A r c ) gauge theories (at zero 
temperature and chemical potential) with content con- 
sisting of Nf copies of massless chiral superfields <&i and 
4>i, i = l,...,Nf, which transform according to repre- 
sentations R and R, respectively, where R denotes the 
representation of the gauge group. We will present vari- 
ous results for arbitrary R and will analyze the following 
specific representations: fundamental, adjoint, and rank- 
2 symmetric and antisymmetric tensor, denoted F, Adj, 
S2, and A 2 , respectively. Using two- and three- loop per- 
turbative results, we study the evolution of the theory 
from the deep UV to the IR and compare these pertur- 
bative results with exact results We investigate 
how the IR zero of the beta function, calculated to a 
certain loop order, and the mass anomalous dimension, 
7 m , of the composite superfield operator product that 
contains the component-field ipij), evaluated at this IR 
zero, calculated to the same loop order, compare with 
exact results. If the theory evolves from the UV to a 
conformally invariant IR phase (a non-Abelian Coulomb 
phase), there is a rigorous upper bound on j m , namely 
7m < 1 dH-pjl. We apply this to a scheme- independent 
calculation of 7,„ and to the two-loop and three-loop val- 
ues of 7 m at the IR zeros of f3, calculated to the same 
order, to calculate perturbative estimates of the minimal 
value of Nf, denoted Nf^ cr , for which the IR behavior 
of the theory is conformal. We carry out this analysis 
in full detail for the case of §i, 3>i in the F + F repre- 
sentation and give briefer analyses for higher representa- 
tions. By comparing these perturbative predictions for 
Nf, C r with exact results, we obtain a measure of the ac- 
curacy of the perturbative analysis. A related way to 
estimate Nf^ cr is from an approximate solution of the 
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Dyson-Schwinger equations for the relevant propagators. 
The comparison of the resultant predictions with exact 
results for the fundamental representation was given in 
the important paper [l6j . Our work complements Ref. 
[l6| . since we analyze the perturbative 7 m rather than 
approximate solutions to Dyson-Schwinger equations. 

There are several motivations for this work. First, it 
is always of fundamental field-theoretic interest to com- 
pare how well a perturbative calculation reproduces an 
exact result. This is especially true in a quantum field 
theory in view of the fact that a pertubative expansion is 
not a Taylor series expansion, with finite radius of con- 
vergence, but instead, only asymptotic, with zero radius 
of convergence [l7|. This is to be contrasted, for exam- 
ple, with high-temperature series expansions in statistical 
mechanics and strong-bare-coupling expansions in lattice 
gauge theory, which are true Taylor series with (at least) 
finite radii of convergence. Secondly, given the great cur- 
rent interest in the nature of the UV to IR evolution of 
asymptotically free gauge theories, as a function of their 
fcrmion content, it is valuable to have a quantitative mea- 
sure of the accuracy of the (semi)perturbative approach 
of calculating the IR zero of the beta function and eval- 
uating the anomalous dimension evaluated at this zero 
of j3. We have previously carried out such a study for 
non-supersymmetric S\J(N C ) gauge theories with various 
fermion contents [18| (see also [19j, whose results are in 
agreement with those in [ijj]). An advantage of mak- 
ing this comparison in a supersymmetric gauge theory is 
that, in contrast to the non-supersymmetric case, one can 
compare the perturbative calculation with exact results 
on the IR properties of the theory. 

This paper is organized as follows. In Sect. II we de- 
fine some notation and recall the relevant coefficients of 
the beta function. Sect. Ill is devoted to a discussion 
of the perturbative expression for 7 m . In Sect. IV we 
review some exact results on the IR phase structure of 
the theory. In Sect. V we give a general discussion of 
perturbative estimates of Nf^ cr . The subsequent sections 
VI- XIII contain our results for the various representa- 
tions considered here. Our conclusions are in Sect. IX. 



II. BETA FUNCTION 
A. General 

The beta function of the theory is denoted /3 = dg/dt, 
where dt = dln/i. In terms of the variable 
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(2-1) 



the beta function can be written equivalently as /3 a 
da/dt = g/3/(2ir), expressed as a series 



da 
~dt 



where £ denotes the number of loops involved in the cal- 
culation of bg and bg = b^j (47r)^. The first two coefficients 
in the expansion (|2.2|) . which are scheme-independent, 
are [H 



3C A 



2T f N f 



and un-m 



b 2 = 6Ci 



4(C A + 2Cf)T f N f 



(2.3) 



(2.4) 



A commonly used regularization scheme for supersym- 
metric theories is dimensional reduction with minimal 
subtraction, denoted DR [3, Hi| (a recent discussion is 
[26j). In this regularization scheme, the coefficient of the 
three-loop term in the beta function is [U 



32] 



-2a h a 1 = -2a ^ b e a 1 , (2.2) 



b 3 = 21C\ + 4(-5C^ - 13CUC/ + AC})T f N } 
+ 4(C A + 6C f )T]N] . (2.5) 

Although the beta function coefficients bi with I > 3 
are scheme-dependent, the use of three- and four-loop 
coefficients in the comparison of the QCD beta function 
with experimental data has shown the value in incorpo- 
rating these higher-loop contributions (33[- The purpose 
in our present work is different from that for QCD; there, 
one wanted to obtain the most precise comparison pos- 
sible with data, in order to extract the value of a s ((i). 
Here, we would like to compare the predictions of a par- 
ticular scheme, namely DR, with the exact results con- 
cerning the infrared behavior of the theory for various 
ranges of Nf, in order to obtain a measure of the accu- 
racy and reliability of the perturbative calculations as a 
guide to the infrared phase structure of the theory. 

Before proceeding, it is appropriate to include several 
cautionary remarks. First, as is well-known, expansions 
such as (|2.2p and (|3.1|) in powers of a are not Taylor 
series, but instead, only asymptotic series, with zero ra- 
dius of convergence. However, a wealth of experience in 
particle physics has shown that if the effective expansion 
parameter (here, (a/ir) times various group invariants) is 
not too large, then the first few terms can provide useful 
information about the physics. Second, the expansions 
for /3 and j m in Eqs. (|2.2p and (|3.1[) are perturbative 
and do not incorporate nonperturbative properties of the 
physics, such as instantons. Instanton effects are absent 
to any order of a perturbative expansion in a, but play 
an important role in a non-Abelian Yang Mills gauge 
theory. Terms arising from instanton effects character- 
istically involve essential zeros of the form exp(— kit /a), 
where k > is a numerical constant. Indeed, instanton 
effects play an important role in the derivation of exact 
results on supersymmetric gauge theories 1 81 . j34| . 

The requirement that the theory be asymptotically free 
means that /3 < 0, which, with the overall minus sign in 
Eq. (|2.2p . is true if and only if £>i > 0. Note that, a 
priori, the condition /3 > could be satisfied with b± = 
if 62 > 0, but this is actually not possible, because the 
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value of Nf that renders b\ = 0, namely [3C 
3C A 



2T 



f 



Nf,, 



(2.6) 



yields a negative value of 62 for any representation i?, 
viz., — YlCApf. Hence, the requirement of asymptotic 
freedom implies 



Nf < N f , max 



(2.7) 



The number Nf, max depends on the representation R, 
and, where necessary for clarity, we shall indicate this 
by writing Nf, max ,R. We shall assume Nf satisfies this 
upper bound in our present study. Although one could 
generalize the analysis to non-asymptotically free theo- 
ries, the coupling a{[i) — > as the energy scale [i — > in 
such theories, so the infrared behavior would be that of 
a free theory. 



B. Zero of the Two-Loop Beta Function 

Since only the two-loop beta function is scheme- 
independent, at the perturbative level, if it does not have 
an IR zero, then, even if such a zero were present at 
the level of three or more loops, it could not be reliably 
considered to be physical. Here we have an even more 
stringent criterion, based on the exact results of Ref. Q , 
which specify, as a function of iV c and the chiral fermion 
content, whether the theory evolves from the UV to a 
conformally invariant IR phase (a non-Abelian Coulomb 
phase) . These results are equivalent to having an ex- 
act beta function and knowing whether it has an exact 
IR fixed point of the renormalization group. If the exact 
analysis does not have an IR zero but the perturbative 
2-loop beta function does have an IR zero, then even 
though the latter is scheme-independent, one would still 
have to reject its prediction, since it differs from the ex- 
act result. In Ref. (3?]], a nonperturbative IR zero of 
the beta function of a non-Abelian gauge theory has also 
been discussed. 

For zero and sufficiently small Nf, the coefficients 62 
and &3 are both positive. As Nf increases, these coeffi- 
cients both decrease. The coefficient 6 2 passes through 
zero and reverses sign from positive to negative at the 
value 



N 



f,b2z 



2T f (C A + 2C f ) 



This value of Nt is less than N 



f 

from the fact that 



Nf, 



(2.8) 



as is clear 



N 



N 



f,b2x 



/,61a 



2Cj_ 

Ca 



< N 



/,61a 



(2.9) 



The two-loop (2£) beta function has a zero away from 
the origin at am,2t = —61/621 i-e-, 



OtIR.2£ 



47T&1 

~6T 



2tt(3C a - 2TfNj) 
2(C A + 2Cf)T f Nf - ZC\ 



(2.10) 



Clearly, for Nf only slightly larger than Nf,b2 Z , aiB.,21 is 
too large for this perturbative result to be trustworthy; 
a necessary condition for it to be reliable is that Nf is 
sufficiently far above Nf,b2z that am,2i is not too large. 
For our analysis below, it will be important whether the 
formal divergence in am^i at 62 = 0, i.e., Nf — Nf,b2z, 
occurs above or below the lower boundary of the IR con- 
formal phase, which is given by Nf, cr in Eq. (|4.5[) . The 
difference is 



N 



f.b2z 



N 



f,cr 



3C A (C A -2C f ) 
AT f {C A + 2Cf) 



(2.11) 



We find that this can be positive or negative. For exam- 
ple, for the fundamental representation, 



Nf, b 2z - N f ,, 



3Nr. 



> for fund. rep. 



2(2JV* - 1) 

(2.12) 

so that 62 = and azR,2i diverges within the IR confor- 
mal phase. In contrast, for the adjoint representation, 

N f}b 2z - N Ur = -- for Adj. rep. , (2.13) 

so that in this case, 62 is nonzero (and negative) all 
throughout the IR conformal phase. For the symmet- 
ric and antisymmetric rank-2 tensor representations, we 
find 



Nf,b2z — Nf, cr 



3N C (N? ± 2N C - 4) 
2(iV c ±2)(3A c 2 ±2A c -4) 

for S 2 , A 2 rep. , (2.14) 



where the upper and lower signs apply for the S2 and 
A 2 representations, respectively. For S2, the numerator 
factor N% + 2N C — 4 vanishes at the unphysical, negative 
value N c = -(1 + v 7 ^ ) and at N c = -1 + V5 ~ 1.236, 
which is less than the minimal non-Abelian value, N c — 
2. Hence, for the 5*2 representation, Nf,b2 Z < Nf, cr for all 
non-Abelian N c and 62 has fixed (negative) sign through- 
out the IR conformal phase. For the A 2 representation, 
N c is restricted to the nontrivial range N c > 3. In this 
A2 case, the numerator factor N^ — 2N C — 4 vanishes at 
N c = 1 + \/h ~ 3.236 (as well as at the negative, un- 
physical value N c — 1 — V5), so that Nf,b2z > Nf, cr for 
the real interval 3 < N c < 1 + y/E, while Nf,b2z < Nf, cr 
for N c > 1 + \/5, i-e., the integer values N c > 4. Note 
that the A 2 representation with A^ c = 3 is equivalent to 
the conjugate fundamental representation. Hence, for all 
representations R for which A2 is distinct from the fun- 
damental, 62 has fixed (negative) sign throughout the IR 
conformal phase. 

Given that Nf < Nf, max to maintain the asymptotic 
freedom of the theory, this am, 2 £ is positive and hence 
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physical if and only if Nf lies in the range Nf^ 2z < Nf < 



3C1 



2T f (C A + 2C f ) 



<N f < 



3C A 
2T f 



(2.15) 



We will thus focus on this interval for Nf . The zero of the 
two-loop beta function at a — am_ 2 i is either an approx- 
imate or exact infrared (IR) fixed point (IRFP) of the 
rcnormalization group. If the gauge interaction sponta- 
neously breaks the global chiral symmetry of the theory 
via the formation of a bilinear matter (chiral) superfield 
condensate, then this IR zero is only approximate, since 
in this case the matter superfield picks up a dynamically 
generated mass £ and as the scale /j, decreases below S, 
one integrates out the matter superfields in defining the 
effective low-energy field theory. Consequently, as the 
theory evolves further into the infrared, the beta function 
becomes that of the pure supersymmetric gauge theory 
without these matter superfields, and hence a(/x) evolves 
away from the approximate IR fixed point. 



C. Zeros of the Three-Loop Beta Function 

To three-loop order, the beta function formally has 
two zeros away from the origin, given by the equation 
bi + b 2 a + b-30 2 = 0, where a was given in Eq. (|2.ip . The 
solutions, in terms of a = Ana, are 



2n 



-b 2 ± x /b 2 2 -4b 1 b 3 



(2.16) 



Only the physical, smaller one of these two solutions will 
be relevant for our analysis, and we label it as am,zi- As 
discussed above, the requirement that the two-loop beta 
function has an IR zero means that Nf is in the interval 
(j2~T5|) where bi > and b 2 < 0. 



III. ANOMALOUS DIMENSION 7m 

The anomalous dimension 7 m describes the scaling 
properties of the quadratic superfield operator product 
containing the bilinear product ip T Cip, or equiva- 
lently, ijjip, of component fermion fields. If one has an 
input mass m for -0, then, with our definition, 7m = 
—dhim/dt, where t = ln/i. Since we are studying the 
evolution from the UV to the IR conformal phase, we 
do not put in such a bare mass m here, since, if we did, 
then as m decreases below m, these fields would be inte- 
grated out as the theory evolved deeper into the infrared 
and the IR behavior would be that of a supersymmetric 
SU(iV c ) theory with just gluons and gluinos. For nota- 
tional simplicity we will often suppress the subscript m. 
This anomalous dimension can be expressed as a series 
in a or equivalently, a: 



7n 



(3.1) 



where eg = ci/{AnY is the £-loop series coefficient. 
We denote 7„£ as the n-loop value of 7m , i.e., j n £ = 

The coefficients C£ have been calculated to three-loop 
order. The one-loop coefficient c\ is scheme- independent: 



(3.2) 



The higher-loop coefficients ci with i > 2 are scheme- 
dependent. In the DR scheme, c 2 and C3 are (28l. l32j 



c 2 = 4Cf(-2C f + 3C A - 2T f N f ) 



(3.3) 



and 



J 



c 3 = 8C f 



AC} + 3C A (C A - C f ) + T f Nf (-8 + 12C(3))C> + (1 - 12C(3))^ 



2T f N f 



(3.4) 



r 



where £(s) is the Riemann zeta function, with £(3) = 
1.20205690.. As Nf approaches Nf :inax from below, b\ — > 
with nonzero b 2 and hence aiR — > 0; since the pertur- 
bative calculation expresses 7 m in a power series in a, it 
follows that 7 m — > as N 



f 



N 



f,max • 



The n-loop value of 7 m at the IR zero of f3, calculated 
to the same loop order (IR fixed point of the renormaliza- 
tion group), is obtained by setting a = a>iR tn i in j n £ (a) 
and is denoted 



where the dependence on the chiral superfield represen- 
tation R is implicit. Thus, at the two-loop level, 



JiR,2i = a(c\ + c 2 a)\ a=aiR: 
h(-cib 2 + c 2 bi) 



(3.6) 



(3.5) Explicitly, 
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liR,: 



C f (3C A - 2T f N f )(2TfN f - C A )(2T f N f - 3C A + 6C» 



[2(C A 



2C f )T f N f - 3Cff 



(3.7) 



r 



Thus, 7/i? ; 2£ has, formally, three zeros and one pole. One 
of the zeros occurs at Nf = Nf tTnax , as given in Eq. (|2.6[) . 
The second zero occurs at 



iV 



Ca 
2T f 



N 



f 7 max 



2N 



/>cr 



(3.8) 



Because this lies below the exact Nf jCr , it is not directly 
relevant for our current analysis. Furthermore, for the 
representations of interest here, it also lies below Nf t b2z> 
and hence is not present where the theory has a two-loop 
zero in /3. The third formal zero in jir^i occurs at 



N f = 



3(C A - 2C f ) 
2T t 



Nf,r. 



3C 



if& • (3-9) 



This third zero occurs for Nf less than Nf cr (and, for 
some representations, at negative Nf), and hence also 
will not be relevant for our analysis. The pole in jib.,21 
occurs at Nf^sz, an( i is a consequence of the pole in 
aiR,2t where 62 = 0. Clearly, the two- loop calculation 
of 7 m ceases to be reliable for Nf values greater than 
Nf.b2z, so this pole is obviously an unphysical artifact. 
Thus, over the range of interest here, 'yiR^e increases 
monotonically above zero as Nf decreases below Nf max . 

In the procedure described above, one evaluates the 
n-loop expression for 7 m at the IR zero of j3, calculated 
to the same n-loop order. For this procedure, one nec- 
essarily uses the two-loop j3 function, since an IR zero 
only appears at this loop level, and also the two-loop or 
higher- loop expressions for j m . Since the coefficients q 
for t > 2 are scheme-dependent, this process necessar- 
ily entails scheme-dependence. We thus also present an 
alternate perturbative estimate for 7 m , which has the ad- 
vantage of preserving scheme-independence but the dis- 
advantage of mixing different orders of perturbation the- 
ory. For this alternative estimate, we use only scheme- 
independent (SI) inputs, and hence evaluate the one-loop 
expression for j m at the two-loop IR zero of /?, obtaining 



lm,si = cxaiR^t = ciarR t 2£ 



2Cf(3C A - 2T f Nf) 
2{C A + 2Cf)TfN f - iC\ 



(3.10) 



IV. REVIEW OF SOME EXACT RESULTS 

Since we will compare our perturbative results with 
certain exact results, a brief review of these is appropri- 
ate. For a vectorlike SU(7V C ) gauge theory with M = 1 



supersymmetry and Nf copies of massless chiral super- 
fields $j and <i>i in the fundamental and conjugate funda- 
mental representation, respectively, exact results on the 
phase structure and corresponding properties of the the- 
ory in the infrared were derived by Seiberg Q. These 
results were subsequently generalized to theories with 
gauge groups SO(7V c ) and Sp(iV c ) in Q and [l(| (reviews 
include (Tl|)- A further gen eralization to arbitrary repre- 
sentations was given in fl2| • In our present work we will 
focus on the comparison of perturbative estimates and 
exact results concerning the minimal value of Nf (for a 
given chiral superfield content), denoted Nf tCr such that, 
for Nt > Nf >cr , the theory evolves from the UV to the 
IR in a chirally symmetric manner, so that the IR theory 
is a conformal, non-Abelian Coulomb phase. This value, 
Nf iCr , is often called the lower end of the conformal phase 
or conformal "window" (with the upper end, Nf^ max , de- 
termined by the requirement of asymptotic freedom). We 
shall carry out this comparison at the two- and three-loop 
level. 

We recall how, for a given R, the conformal region in 
Nt is determined. A crucial tool in determining Nf jCr , 
the lower end, as a function of Nf , of the IR conformal 
phase, is the existence of an exact relation between the 
beta function and the mass anomalous dimension, j m . 
This relation is embodied in the following form for the 
beta function of the theory with a vectorlike massless 
chiral superfield content consisting of Nf cop ies of the 
representations R + R of the gauge group [34 l3l| : 



da a 2 
^ a = l± = ~2^ 



b 1 -2T f Nf lm (a) 



1 



C A a 
2ir 



The IR zero of j3 a is determined by the condition 
3C A — ZTfNf _ Nf t7nax 



7n 



2T f Nf 



Nr 



1 



(4.1) 



(4.2) 



Let us assume that the theory flows to an exact IR fixed 
point, i.e., that Nf is in the sub-interval of (|2.15l) in 
which, as the theory evolves down from the UV to the 
IR, no spontaneous chiral symmetry breaking takes place. 
Given that the theory has evolved down to an (exact) in- 
frared fixed point, the resultant theory at this IRFP is 
conformally invariant. 

It is a special property of a conformally invariant field 
theory (whether supersymmetric or not) that the full di- 
mension of a spinless operator (other than the identity) 
must be larger than unity in order that the theory not 
contain a ny neg ative-norm states, which would violate 
unitarity |13l4l5l |. Specifically, for the dimension D m of 
the bilinear operator <&i<&i (with no sum on i) for any 
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i = 1, ...,Nf in the present theory, this is the inequality 

D m > 1 . (4.3) 

In terms of its component scalar and fermion fields 4>i 
and ipi, the chiral superfield is expressed as <£>i = 4>i + 
8ipi + 88 Fi , where 9 is a Grassmann variable and Ei is 
an auxiliary field. Thus (for any i), the term <&i$j yields, 
as the (holomorphic) term bilinear in component fermion 
fields, 98ipi4>i- Taking into account that the dimension 
of 8 is —1/2, the free-field dimension of ipi4>i (for any 
i) is 3, and using our definition of "f rn , it follows that 
D m = 2 — 7 m , so that the bound (14.31) is equivalent to 
the following upper bound on j m : 



In 



< 1 



(4.4) 



This may be contrasted with the situation in a non- 
supersymmetric SU(-/V C ) theory. There, the bound that 
the full operator dimension of ipiipi be larger than 1 im- 
plies that 7 m < 2, as we noted in Eq. (4 2) of (l8j (equiv- 
alent to the bound from Eq. (4.1) of [HI]). The more 
stringent upper bound (|4.4I) on j m in the supersymmet- 
ric theory is due to the fact that the fermion field ipi is 
part of a chiral superfield and the holomorphic fermion 
bilinear resulting from the quadratic <J?j<&j product car- 
ries with it a 98 factor. 

We next assume that, in the relevant range of Nf where 
the theory evolves from the UV to an IR-conformal phase, 
7 m evaluated at the IR fixed point, a = am, increases 
monotonically as Nf decreases below Nf^ max . This as- 
sumption is satisfied by 7 m as calculated in a scheme- 
independent manner, as will be discussed further below. 
The inequality (14.41) then implies that Nf tCr , the value 
of Nf below which the theory cannot be conformally in- 
variant, is bounded below as Nf tCr > 3C A /(4Tf). The 
application of duality relations provides strong evidence 
that this inequality is saturated @, [l2j and hence that 



Nf,, 



2>C A 



(4.5) 



We refer to this as an exact result, although, as we have 
indicated, there are some nonrigorous steps in its deriva- 
tion. Note that 



N f,cr = 



N 



f,max 



(4.6) 



Thus, the theory evolves from the UV to an IR fixed 
point in the conformal phase if and only if iV/ jCr lies in 
the interval 



3CA <N f < 3CA 



2T 



I 



(4.7) 



(The marginal value Nf = 2>C A /{ATf) itself is not in this 
conformal phase @,H|.) For both Eqs. (|4~5j) and (j4~7] . 
it is understood that, physically, Nf must be an integer 
36]. Thus, the actual values of Nf in the conformal 
phase are understood to be the integers that satisfy the 
inequality (|4.7|) . 



PERTURB ATIVE ESTIMATES OF N f . cr FOR 
GENERAL R 



As discussed above, although a perturbative calcula- 
tion is not exact, one gains valuable information by com- 
paring it with exact results. We carry out this com- 
parison here for a general representation R, using per- 
turbative estimates for Nf^ cr , the lower boundary of the 
IR conformal phase. For this purpose, we utilize "fm,si 
and ^im,2i- With a monotonic increase in 7 m as Nf de- 
creases below Nf^ max , we can then calculate a perturba- 
tive estimate for Nf^ cr by assuming that 7 m saturates the 
inequality l|4.4p as Nf decreases through Nf^ cr . (Here, 
again, we are implicitly analytically continuing Nf from 
physical integer values to real numbers.) Setting the per- 
turbative 7,„ = 1 and solving for the value of Nf at which 
this happens yields the corresponding perturbative esti- 
mate of Nf jCr . Since a perturbatively calculated expres- 
sion for 7 m is not, in general, equal to the exact j m , one 
does not expect these estimates to agree precisely with 
the exactly known values for Nf^ cr for the various repre- 
sentations. However, this comparison gives quantitative 
insight as to the accuracy of the perturbative calcula- 
tions. 



A. Estimate Using 7/h,s/ 

The scheme-independent perturbative result for j m , 
Iir.si, increases monotonically as Nf decreases from 
its maximal value (|2.6p and reaches the rigorous upper 
bound as Nf decreases through the value 



N 



f,cr,SI 



3C A (C A + 2C f ) 
2T f (C A +ACf) 



(5.1) 



This is larger than the exact value of Nf_ cr , as is evident 
from the difference 



N 



f,cr,SJ — JV /,cr 



Nf. nr = 



4T f {C A +ACf) 



> (5.2) 



or the ratio 



N 



f,cr,SI 



Nf,, 



= 2 



C A + 2C 



C A + 4C f 



> 1 



(5.3) 



This difference between the scheme-independent pertur- 
bative estimate of the lower boundary of the conformal 
phase, Nf^ cr .si, and the exact lower boundary, Nf tCr , 
provides one quantitative measure of the accuracy of 
perturbation theory. Our conclusion from this compar- 
ison is that perturbation theory slightly overestimates 
the value of this lower boundary and hence underesti- 
mates the size of the conformal phase as a function of 
Nf. Related to this, as Nf decreases below Nf_ cr: si 
toward the exact lower boundary of the IR conformal 
phase at Nf^ cr , "/ir^si continues to increase. In this 
regime, its behavior is unphysical since it violates the 
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rigorous bound (|4.4[) . This happens for both represen- 
tations where Nf^z > Nf cr and representations where 
Nf, b 2z < N fyCr . Formally, 



1IR,S1 



2C f 



2Cf — Ca 



at N f = N 



(5.4) 



B. Estimate Using yin t 2£ 

Setting the two-loop result for the anomalous dimen- 
sion, jiR^e, equal to the rigorous upper bound, unity, 
we derive the corresponding two-loop perturbative pre- 
diction for Nf jCr . The equation jiR^e = 1 is a cubic 
equation in Nf, which yields the resultant estimate for 
Nf^ cr , together with two other roots that are not of direct 
relevance. Formally, 



llR,2l 



C f (AC f -C A ) 
2(2C> - C A ) 2 



at N f = N f>cr 



(5.5) 



but as with "fm sii this is only formal, since this pertur- 
bative result generically violates the upper bound (|4.4p . 
We comment below on the situation at the three-loop 
level. We proceed to present results for the various rep- 
resentations of interest here. 



VI. CHIRAL SUPERFIELDS IN THE 
FUNDAMENTAL REPRESENTATION 



IR Zeros of the Beta Function 



Two-Loop Ana, 



In this section we consider the case where the theory 
has Nf copies of massless chiral superhelds $j and $i, 
i = 1, ...,Nf, transforming according to 



4\ 



*i: F, i = l,...,N f , 



(6.1) 



i.e., the fundamental plus conjugate fundamental, repre- 
sentation of the gauge group. The requirement of asymp- 
totic freedom implies 

N f < 3N C . (6.2) 

For this case, the exact result (|4.5p on the value of Nf at 
the lower boundary of the conformal phase in the infrared 
is 



3N C 
2 



(6.3) 



where it is understood that this is only formal if N c is 
odd, since Nf tCr must be an integer. Thus, the IR con- 
formal phase is given, from Eq. (|4.7j) . as 



3iV e 
2 



< Nf < 3N C 



(6.4) 



TABLE I: Values of Nf : t,i z = Nf )ma x> N f,b2z, and Nf : t, 3z for the 
supersymmetric SU(A r c ) theory with Nf chiral superfields <t>i, <t>i 
in the F and F representations, respectively. We also list the exact 
value of Nj cr . These results are given for the illustrtive values 
2 < N c < 5. 



N c 


N f,cr 




N f,b3z 




2 


3 


3.43 


3.09 


6 


3 


4.5 


4.76 


4.27 


9 


4 


6 


6.19 


5.55 


12 


5 


7.5 


7.65 


6.85 


15 



Physically, Nf must be a (non-negative) integer, so the 
actual physical values of Nf in the IR conformal phase 
for 2 < N c < 5 are N c = 2 : N f = 4, 5; N c = 3 : 



N 



f 



5, 6, 7, 8; N c = 4 : N f = 7, 8, 9, 10, 11;, and 
N c = 5 : Nf — 8, 9, 10, 11, 12, 13, 14. 

Evaluating Eq. (I2.8[) . we find that 62 reverses sign from 
positive to negative as Nf increases through the value [36[ 



N 



3Nr. 



/,62a 



(6.5) 



The interval of Nf values in Eq. (|2.15j) where the two- 
loop beta function has an IR zero is therefore |36| 



3N C 



2-N c 



<N f < 3N C 



(6.6) 



Numerical values of Nf :Cr , Nj t b2z, Nf t b3z, and Nf^u 



N 



f,max 



are listed in Table Q] for the illustrative values 
2 < N c < 5 . As discussed before in connection with 
Eq. (|2.12p . the value of Nf, cr in Eq. (|6.5p is greater 
(for all N c ) than the exactly known lower boundary of 
the conformal phase in Eq. (|6.3j) . As iV c — > 00, Nf^z 
asymptotically approaches (Z/2)N C from above. 

Since Nf^2z > Nf_ cr , it follows that the two-loop beta 
function only has a (perturbative) infrared zero for Nf 
values in the interval where the theory is conformally in- 
variant. This is different from the non-supersymmetric 
SU(A r c ) gauge theory, in which the two-loop beta func- 
tion may have an IR zero for values of Nf less than the 
estimate, from the Dyson-Schwinger equation, of Nf jCr , 
i.e., in the phase where the theory has spontaneous chiral 
symmetry breaking. (Because of this SxSB, this IR zero 
is only approximate.) 

For our supersymmetric theory with the F + F chiral 
superfield content of Eq. (|6.ip . the general formula (|2.10p 
for the IR zero of the two-loop beta function reduces to 



2tt(37V c - Nf) 



(2N C - Nc )Nf - 3AT? 



(6.7) 



This decreases monotonically from arbitrarily large val- 
ues (where, of course, the perturbative beta function does 
not apply reliably) to zero as Nf increases throughout 
the interval (I6.6p . Numerical values of a-m^i are listed 
in Table [TTI for the illustrative cases 2 < N c < 5. 
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TABLE II: Values of the IR zero of the beta function in the super- 
symmetric S\J(N C ) gauge theory with Nf pairs of chiral superfields 
in <E>i, <&i in the fundamental and conjugate fundamental repre- 
sention, respectively, calculated at n-loop order, and denoted as 
a lR nl- Results are given for the illustrative values 2 < JV C < 5. 
For each N c , we only give results for the integral Nf values in 
the interval 1 12, 151 1 where the theory is asymptotically free and the 
two-loop beta function has an infrared zero. 



AT 


AT 


ClIB.,21 




n 
z 


4 


O.ZO 


Z.OO 


9 





11/1 
1.14 


u.oyo 


n 


r 




lo or 

lo.oo 


o.Uo 


Q 

o 







1 /in 
1.4U 


Q 


*7 

( 


n ono 

u.yyz 


U. / o4 


Q 





r\ q /i q 


U.oUo 


4 


7 


5.03 


1.64 


4 


8 


1.795 


0.984 


4 


9 


0.867 


0.615 


4 


10 


0.426 


0.357 


4 


11 


0.169 


0.158 


5 


8 


12.94 


1.90 


5 


9 


2.86 


1.13 


5 


10 


1.37 


0.765 


5 


11 


0.766 


0.528 


5 


12 


0.442 


0.353 


5 


13 


0.240 


0.212 


5 


14 


0.101 


0.0963 



It is often of interest to consider the 't Hooft limit 
N c — > oo with g 2 N c fixed and finite. For the present 
F + F superfield content it is also natural to consider 
taking Nt — > oo with the ratio 



N 



f 



Nr. 



(6.8) 



fixed and finite (sometimes called the Veneziano limit). 
In this limit, the relevant interval for r where the two- 
loop beta function has an IR zero is thus 

3 
2 

Here. 

2tt(3 - r) 



< r < 3 



(6.9) 



OtlR2iN c 



(6.10) 



2r-3 

which decreases monotonically to as r increases through 
the interval 3/2 < r < 3. 

2. Three-Loop Analysis 



For the present case, Eq. (|6.1|) . the general result in 
Eq. (12.51) for the three-loop coefficient 63 takes the form 

63 = 217V C 3 + (9 - 21N 2 C + 2N~ 2 )N f 



3.11) 



For small Nf, 63 is positive. As Nf increases, 63 passes 
through zero and reverses sign from positive to negative. 
To investigate this, one solves the equation 63 = for Nf. 
Since 63 is a quadratic function of Nf , there are formally 
two solutions to this equation, namely 



N f , b3 z,± = [2N C (AN 2 - 3)]- 1 - 9N 2 - 2 



3.12) 



For N c — 2, this gives Nf t b3 Z ,- — 3.09 (to the indicated 
accuracy), slightly above (3/2)^ = 3, while for N c > 3, 
we find that Nf^3 Z < (3/2)N c . As examples, for A^ c = 3, 
N f , b3z ,- = 4.27 < (3/2)N c = 4.5, while for N c = 4, 
Nf, b 3z - = 5.55 < (3/2)N c = 6, and so forth for higher 
values of N c . For large N c , 



(3/2)N c 



1 

12 



21 - 



O 



1 

N? 



3.13) 



The numerical value of the first term is approximately 
0.896. As Nf increases past the larger value Nf t b3 Z ,+, 
&3 vanishes and reverses sign again, becoming positive. 
However, this larger zero is not relevant to our analysis, 
since for all N c > 2, Nf,b3z,+ > 3iV c , i.e., this occurs for 
Nf beyond the upper limit imposed by the constraint of 
asymptotic freedom. For example, for N c — 2, Nf^ z ,+ — 
8.38, which is greater than Nf^ max — 6, and so forth for 
larger values of A^ c . For large N c , we have 



3AL 



1 

24 



O 



1 



3.14) 



The numerical value of the first term is 1.30. 

Since we restrict here to the interval (|2.15j) where the 
two-loop beta function has an IR zero, we observe that 
for physical, integral values of Nf in this interval, 63 is 
always negative. Note that for N c = 2, 63 vanishes and 
changes sign from positive to negative as Nf increases 
through Nf — 3.09, but this value of Nf is less than the 
value Nf^2z — 3.43 where the beta function first has an 
IR zero for this N c . Hence, for all N c > 2 and for integral 
values of Nf in the interval (|2.15[) where there is an IR 
zero at the two-loop level, 63 < 0. 

We list values of a IRt3l for 2 < N c < 5 in Table HH 
together with the values of am^i already given. As is 
evident from this table, 



(6.15) 



This is the same trend that we found in |Ts] for a 
non-supersymmetric SU(A r c ) with Nf copies of massless 
fermions in the fundamental representation. 
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In the large-A^, large-A/ limit with Nt = rN c 
calculate 



N c . = 



4tt 3 - 2r + V4r 3 - 29r 2 + 72r - 54 



21(r-l)-4r 2 

(6.16) 

The right-hand side of Eq. (|6.16|) decreases monotoni- 
cally from 4tt ~ 12.57 to as r increases from 3/2 to 3. 
Note that the denominator in Eq. (|6.16p . 21 (r— 1) — 4r 2 , 
is positive-definite in the relevant interval 3/2 < r < 3; 
it has zeros at r = (1/8)(21 ± \/105), i.e., approximately 
3.906 and 1.344. In the numerator, the first part, 3 — 2r, 
is negative-definite in this interval, but is smaller than 
the square root. 



J 



B. Values of j m at IR Zero of j3 

1. Coefficients 

Evaluating Eqs. (|3.2[) - (|3.4j) for our case (|6.1|) . we have 
2(A 2 - 1) 



Nr. 



3.17) 



c 2 = 2{N 2 - l)N- 2 {2Nj + 1 - N c N f ) , (6.18) 



and 



c 3 = 2{N 2 C - 1)N~ 



5iV c 4 - N 2 + 2 + N C (-3N 2 + A)N f - 6C(3)A C (A 2 + l)N f - N 2 c N 2 f 



V 



(6.19) 



where £(s) is the Riemann zeta function. As Nf ap- 
proaches Nf <max from below, b\ — > with nonzero 6 2 
and hence am — > 0; since the perturbative calculation 
expresses j m in a power series in a, it follows that as 
7 m -> as iV/ -> N f . max . 

To get an analytic understanding of the behavior of 7,„, 
we study the signs of the coefficients C£ with £ = 1, 2, 3. 
Since the one-loop coefficient c\ is positive, j m increases 
from zero as Nf decreases just below Nf^ max . In the 
conformal phase the two-loop coefficient ci may be either 
positive or negative, depending on Nf. This coefficient 
C2 vanishes at 



N f = 2N C + N- 1 =N 



f,c2z 



(6.20) 



and 



(-2 



> for N f < N 



*f,c2z 



c 2 < for N fiC2z < Nf < TV 



f,max 



3N C .(6.21) 



The three-loop coefficient C3 is a quadratic function of 
Nf and vanishes at two values of Nf = 0, namely 



N 



f,c3z,fund,± 



4 - 37V 2 - 6(7V 2 + 1)C(3) ± VR^~z 



2N r 



where 

R c3z = 297V C 4 - 28N 2 + 24 + 12(7V C 2 
+ [6(N 2 + 1)] 2 C(3) 2 . 



3.22) 



1)(37V 2 -4)C(3) 



3.23) 



For N c = 2, iV f,c3zjund,- is equal to —22.88 and hence is 
unphysical, while Nf >c s z Jund,+ = 0.8522, so that in the 
physical range, C3 is negative for all positive integral val- 
ues of Nf. For N c = 3, Nf, c3z j un d,- is equal to —33.05 



and hence is again unphysical, while Nf tC 3 Z ,fund,+ = 
1.338, so that in the physical range, C3 is positive for 
Nf = 1, but negative for Nf > 2, including all of the 
interval of interest here. This qualitative behavior con- 
tinues to hold for higher values of N c , as is evident from 
the Taylor series expansion 



N 



f,c3z,fund,- 
~Nr. 



+ ^36C(3)(C(3) + 1) 



-3(1 + 2C(3)) 



29 



3.24) 



The numerical value of the constant term is 0.468 to the 
indicated accuracy. Thus, as N c — !• 00, Nf iC a z jund,+ ~ 
0.468^, which is less than the exact value of the lower 
boundary of the conformal phase, Nf :Cr = (3/2)JV c , given 
in Eq. (|6.3|) . Thus, a general characterization of the 
three-loop 7 TO is that c\ is positive for any representation 
R, and for R = F; c 2 is positive in the lower part of 
the conformal phase but negative in the upper part, as 
specified by Eq. (|6.21[) . and C3 is negative throughout 
all of the Nf interval of interest, including the conformal 
phase. 



2- jir,si 

For the chiral fermion content (|6.ip . we calculate 



JlR,SI 



(N 2 -1)(3N C -N f ) 
(2iV 2 - l)N f - 3iV| 



3.25) 



As Nf decreases below its maximal value, 3iV c , "fm t si 
increases monotonically. We list values of Jir^i hi Table 
IIIII for the illustrative values 2 < N c < 5. For each N c , 



10 



TABLE III: Values of the anomalous dimension jm in the SXJ(N C ) 
supersymmetric gauge theory with Nf copies of massless chiral 
superfields <J>i, &i in the F and F representations, calculated to the 
n-loop order in perturbation theory and evaluated at the IR zero of 
the beta function calculated to this order, (Xir n i> f° r € = 2, 3. We 
denote these as ~fiR,„i = 7ni( a lR.,ne)- Results are given for the 
illustrative values 2 < N c < 5. For sufficiently small Nf > Nf^2 z 
for each N c , Q-jr 21 is so large that the formal values of ^jr 21 
and /or ■yjR 31 are either larger than unity or negative and hence are 
unphysical. We indicate this by placing these values in parentheses. 



N c 


Nf 


JIR^I 


liR.21 


1IR,U 


2 


4 


(1.500) 


(1.875) 


(-1.68) 


2 


5 


0.273 


0.260 


0.0802 


3 


6 


(1.14) 


(1.22) 


(-0.730) 


3 


7 


0.421 


0.399 


0.0584 


3 


8 


0.145 


0.139 


0.104 


4 


8 


(1.07) 


(1.11) 


(-0.546) 


4 


9 


0.517 


0.490 


0.0219 


4 


10 


0.254 


0.239 


0.127 


4 


11 


0.101 


0.0970 


0.0835 


5 


10 


(1.04) 


(1.07) 


(-0.475) 


5 


11 


0.585 


0.557 


(-0.0135) 


5 


12 


0.338 


0.317 


0.120 


5 


13 


0.183 


0.173 


0.121 


5 


14 


0.0772 


0.0748 


0.0680 



we omit values in the lower range of Nf that strongly 
violate the bound (|4.4[) . This violation is due to both the 
inexactness of the perturbative calculation of 7 m and the 
fact that the two-loop IR zero of the beta function, am t 2i 
gets arbitrarily large as Nf decreases toward Nf^2 Z . In 
Figs. [l][3]we show plots of Jir.si as a function of Nf for 
the approximate respective subintervals of the conformal 
phase where 7/_r,sj satisfies the upper bound (14. 4ft . Wc 
will discuss below the other curves on these plots. 

Specializing Eq. (15.11) to the case of the fundamental 
representation, we obtain the scheme-independent per- 
turbative estimate of the lower boundary of the confor- 
mal phase, 



N 



f,cr,SI 



3N C (2N? - 1) 
3N? - 2 



.26) 



As noted before, this estimate lies above the actual exact 
lower boundary, which, in the present case, occurs at 
N Lcr = (3/2)JV c : 



N fiCr ,si _ 2{2N 2 C - 1) 



Nf,, 



3Ar 2 _ 2 



3.27) 



This ratio decreases from the value 7/5 = 1.40 at N c = 2 
to 4/3 as N c — > 00 and has the Taylor series expansion 



Nf. cr ,SI _ 4 

~ 3 



N 



f,cr 



o 



1 



as N r 



00 



3.28) 




FIG. 1: Plot of the n-loop fcrmion anomalous dimension, -f„£, 
evaluated at the respective n-loop value of the IR zero of 0, a_rij,n£, 
and denoted as jiRnt, for two-loop and three-loop order, in the 
case of Nf chiral superfields <E>i, <J>i in the F and F representation 
of SU(iV c ) for N c = 2. We use the generic label 7 for the vertical 
axis. At the lower end of the plot, from top to bottom, the curves 
are for (i) fiR 2I7 (") llR Sly an£ i ( m ) lIR 3f ■ The curves involve 
an implicit analytic continuation of Nf from integer values to real 
values; of course, only the integer values are physical. We only show 
the region in Nf where JiR2£ an d 7ir si approximately satisfy the 
upper bound 14.41 . 




FIG. 2: Same as Fig. [T]for N c = 3. 




FIG. 3: Same as Fig. [T]for N c = 4. 
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We next proceed to two-loop and three-loop analyses. 
As we have remarked above, these have the advantage 
of using the same n-loop orders in calculating 7 m and 
oiiR,nii but are subject to the standard caution that they 
involve scheme-dependence. 



3. Two-Loop Ana, 



Evaluating j m calculated to two-loop order at the zero 
of the beta function calculated to the same order for R = 
F, we obtain 



JlR,2£ 



(A 2 - 1)(3JV C - N f )(N f - N c )(N c N f - 3) 



2(-3iVf + 2N*N f - N f f 



(6.29) 



We list values of jir^z in Table IIIII and plot curves of 
"tiR,2i as a function of Nf (analytically continued from 
integer to real values) in Figs. [TJO One sees from these 
results that over the range where the calculations are 
reliable, jiR t 2t is rather close to Jir.si- The anoma- 
lous dimension "fiR^i increases monotonically as Nf de- 
creases from Nf, max = 3N C in the interval (|2.15p . As 
Nf decreases toward Nf^z (given in Eq. (|6.5I) . am^i 
gets arbitrarily large, and this perturbative expression 
obviously ceases to be reliable. As Nf increases to- 
ward Nf, max j und = 3A C , aiR^l -> 0, and hence also 
7/fl,2£ -> 0. 

The condition that jiR^e = 1 is a cubic equation in 
Nf, namely 

A C (A 2 - l)Nf + (4A C 4 - 7A 2 + 5)Nf 

+ 3N C (-7N* + 7A 2 - A)N f + 9A 2 (2A C 4 - A 2 + 1) = . 

(6.30) 

As Nf decreases from Nf^ rnax j unc i = 3N C toward 
Nf,b2zjund, 7m exceeds the rigorous upper bound 7 m < 1 
at a value of Nf which is given as the relevant one among 
the three roots of Eq. (|6.30p . For example, for N c = 2, 
this equation has the physical root Nf = 4.242, together 
with two other roots which are not relevant, namely 
Nf = 2.929 and Nf = -14.00 (to four significant fig- 
ures). The first of these other roots is irrelevant since it 
is below the lower boundary of the IR conformal phase, 
Nf cr given in Eq. (|6.3I) . and the second is irrelevant 
since it is negative and hence unphysical. Thus, insofar 
as one can compare a perturbative 2-loop calculation of 
7 m with an upper bound on the exact j m , one finds that 
this 2-loop prediction for the lower end of the conformal 
phase in Nf for N c = 2 is Nf = 4.242. The ratio of 
this to the exact result Nf^ cr , which is 3 for N c = 2, is 
approximately 1.414. Similarly, for N c = 3, the physi- 
cal root of Eq. ([OUf is Nf = 6.150. (The other two 
roots are 3.983, and —21.22, which are again irrelevant). 
The ratio of this to the (formal, half-integral) exact result 



Nf iCr = 9/2 is 1.367. As N c —> oo, we find that this ratio 
of the physical root of the cubic for Nf jCr divided by the 
exact value Nf >cr = (3/2) A c approaches the same value, 
4/3, as was true of the ratio of Nf jCr: si divided by this 
exact value, as given in Eqs. (|6.27l) and (|6.28[) above: 



Nf,cr ~ 3 + \N? 



3.31) 



We thus find that a 2-loop perturbative analysis of (3 and 
7 m overestimates the value of Nf_ cr somewhat and hence 
underestimates the size of the conformal phase in Nf for 
this case (|6.1I) . This is qualitatively the same as we found 
for the estimate of Nf tCr using Jir^si- 

It is interesting to compare these estimates for Nf^ cr 
from the scheme-independent "fiR t si and the two-loop 
liR,2t with the results obtained via the different method 
of equating cxm^e with the critical value a cr for dynami- 
cal mass generation calculated from the Dyson- Schwinger 
equation for the fermion propagator in Ref. 16]. This 
mass generation, associated with fermion condensation, 
was found to occur at UM 



Nf^cr.DS 



3A^ C (3A C 2 



AN} 



3.32) 



The ratio of this to the exact value for Nf lCr in Eq. (|6.3 



Nf^os _ 3/l-(2/3)JV, 



N fiCr 2Vl-(3/4)JV, 



3.33) 



This ratio is equal to 1.54 and 1.515 for N c = 2 and 
N c = 3 and decreases monotonically to 3/2 as N c — > oo. 
Thus, for a given N c , our estimates of Nf iCr obtained 
from equating 'Jir.si and jiR t 2i to the upper bound of 
unity are slightly closer to the exact result (|6.3p than 
the estimate from the analysis of the Dyson- Schwinger 
equation, and all agree qualitatively, i.e., all are slight 
overestimates of Nf^ cr . 

In the large- N c , large- Nf limit, with Nf — rN c , jiR^e 
has the Taylor series expansion 



liRat 



r(r — 1)(3 — r) 
2(3 -2r) 2 



3(r _l )( 3_ r) 2 / ! 
2(3-2r) 3 A 2 \N* 



(6.34) 



In this limit, the value of r, and hence Nf, where r )iR2i 
I is given as the relevant one among the three roots of 
the cubic equation (r — 2)(r 2 — 6r — 9) = 0, namely r = 2. 
The other two roots are — 3 + 3-\/2 and — 3 — 3\/2; the first 
of these has the value 1.243 and hence is below the value 
r = 3/2 (cf. Eq. (|6.5[) ) where, for increasing r, an IR zero 
of the beta function first appears (which coincides with 
the exact value of the lower boundary of the IR conformal 
phase in this limit), and the second is negative and hence 
obviously unphysical. Thus, in this large- N c limit, with 
r = Nf/N c fixed and finite, the 2-loop analysis predicts 
that Nf <cr = 2N C , which is 4/3 times the exact result of 
Eq. (I6.3[) . in agreement with our analysis in Eq. (I6.31[) . 
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4- Three-Loop Analysis 

In the same manner, we evaluate the three-loop ex- 
pression for 7 m at the 3-loop value of the IR zero of the 
beta function am^i- Since the analytic formulas are 
somewhat complicated, we will restrict ourselves to giv- 
ing numerical results for N c = 2 through N c — 4 in Table 
Mil and Figs. Q][3l In contrast to Jir^si an d 1ir,21, we 
find that ^ir^i does not increase monotonically as Nf 
decreases below Nf max . Instead, it reaches a maximum 
well below unity in the interior of the conformal phase 
and then decreases, vanishing and becoming negative. 
Because of this behavior, we cannot use our procedure 
of setting the perturbative expression for 7 m equal to 1 
and then solving for Nf for ^ir^i- Since this behavior 
of 1ir,31 clearly differs from the behavior of the scheme- 
independent jir^si, it may reflect scheme-dependence. 
It also shows that as Nf decreases toward the lower end 
of the IR conformal phase and aiR increases to larger 
values, a perturbative calculation becomes less reliable. 
In general, a necessary condition for these perturbative 
calculations to be reliable is that inclusion of the next 
higher-loop order term should not drastically change the 
qualitative behavior. Imposing this condition for the 
comparison of jir^i and ^yiR^e, we may obtain an es- 
timate of the interval in Nf where the calculation could 
be reasonably reliable. For N c = 2, we find that this 
interval plausibly includes Nf = 5 but does not include 
N f = 4. For N c = 3, this interval includes Nf = 8 
but not lower values of Nf. In view of this behavior of 
liR,3t, one must view the three-loop results with appro- 
priate caution, recognizing that perturbative calculations 
become less reliable as the coupling becomes stronger. 



VII. SUPERFIELDS IN THE ADJOINT 
REPRESENTATION 

A. Beta Function 

The adjoint representation is self-conjugate, so here a 
theory with Nf copies of a massless chiral supcrficld con- 
tent consisting of $i and <&i, i — l,...,Nf, is equivalent 
to a theory with Nf = 2Nf copies of We shall thus 
consider half-integral values of Nf as physical here. The 
beta function coefficients are 

h = N c (3 - 2N f ) , (7.1) 

b 2 = -%Nt{2Nf - 1) (7.2) 

and, in the DR scheme, 

& 3 = -7iV c 3 (2JV)-l)(3-2iV» . (7.3) 

Note that b 3 vanishes at the same (formal, non-integral) 
value of Nf at which b\ vanishes, namely Nf = 3/2. The 



condition that the theory be asymptotically free, i.e., that 
&i > 0, implies the upper bound Nf < 3/2. 



N f < g = N f,max 



(7.4) 



For the &i content, this only allows the choice Nf = 1, 
while for the reduced content consisting only of $i, this 
allows N' f = 1 and N' f = 2. (Note that the N' f = 1 
theory has been solved exactly [38j| . ) 

From Eq. (|4.5j) . the lower boundary of the IR confor- 
mal phase is given formally by 



N 



f.cr 



(7.5) 



or equivalently, N'f cr = 3/2. Since neither of these is an 
integer, they must be regarded only as quantities defined 
via a requisite analytic continuation of the theory in Nf 
or Nf away from the integers to the real numbers. With 
this understanding, the IR conformal phase is thus given 
by 



3^3 
- < N f < - 
4 1 2 



(7.6) 



Hence, with the &i superfield content, the only inte- 
ger value of Nf allowed by the requirement of asymptotic 
freedom, namely Nf = 1 yields an IR conformal phase. 
For the theory with just the $i superfield, Eq. (|7.6[) reads 
3/2 < N'f < 3, so for N' f = 2 (N' f = 1) the theory evolves 
into the infrared in a conformal (nonconformal) manner, 
respectively. 

In the theory with $, $ superfield content, the two- 
loop /3 function coefficient bi is negative for the only 
relevant value of Nf, namely Nf = 1. In the reduced 
theory, b 2 = -6N*(N f - 1), which is zero for N' f = 1 
and negative for N'f — 2. Thus, at the two- loop level, 
the IR zero of the $ functions occurs at 



&IR,2£ 



2tt(3 - 2N f ) 
3N c (2N f - 1) 



(7.7) 



At the three-loop level, j3 has two zeros away from the 
origin, at 



a 

47T 



-3± 



2(14Ar|-33JV / +27) 



7iV c (3 - 2N } ) 



(7- 



The — sign choice yields an unphysical, negative result, 
so the physical three-loop IR zero of the beta function is 
given by Eq. (|7.8|) with the + sign. We denote this as 
o-ir,u (suppressing the Adj for simplicity). This uir,m 
exhibits unphysical behavior, vanishing at Nf =9/8 and 
becoming negative in the range 9/8 < Nf < Nf, max — 
3/2. One could take the point of view that this precludes 
a reliable three-loop perturbative analysis of this case. 
However, we will at least give results for the one case for 
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which the theory has an IR fixed point, namely Nf = 1. 
For Nf — 1, we have 



OiIR,2i,N t =l 



and 



&iR,3e,N f =i 
so that, for this value of Nf, 

OtIR.M 

0-IR,2(. 



2tt 

= ZN C 

47T 

7N~ C 



(7.9) 



(7.10) 



(7.11) 



independent of N c . This is the same trend that we found 
for the case of matter superfields in the F + F repre- 
sentation, i.e., the value of the IR fixed point calculated 
to three-loop order is somewhat smaller than the value 
calculated to two-loop order. 



B. Anomalous Dimension 

For this theory with chiral superfields in the adjoint 
representation, the coefficients in Eq. (|3.ip are 



C-2 



ci = 4N C 



-4N?(2N f - 1) 



and 



C3 



-8N?(N f + 4)(2JV / -1) 



(7.12) 



(7.13) 



(7.14) 



From our general result (I3.10|) . we calculate the scheme- 
independent anomalous dimension 



Jir,SI — 



2(3 - 2N f ) 
3(2N f - 1) 



(7.15) 



This increases monotonically as Nf decreases from its 
maximal to its minimal value in the IR conformal phase. 
1ir,SI increases through its upper limit of unity as Nf 
decreases through the value Nf >cr> si- Evaluating our gen- 
eral formula in Eq. (|5.ip for the present case of the ad- 
joint representation, we obtain 



Nf,er,SI 



9 

10 



(7.16) 



As is true for general R, this is larger than the exact 
result, which in the present case is Nf tCr = 3/4. (This 
exact value is only formal, since it is non- integral.) The 
ratio (|5.3[) here is 6/5 = 1.2. As Nf decreases from 9/10 
to Nf cr — 3/4, Jir^si increases from 1 to 2, exhibiting 
unphysical behavior. 

Evaluating Eq. (|3.6I) for the present case of the adjoint 
representation, we find 



lIR.U 



(3-27V / )(2A / + 3) 
9(2N f - 1) 



(7.17) 



Setting this equal to 1, we derive another perturbative 
estimate of Nf cr , which is the positive root of a quadratic 
equation, 

Nf = H-3 + VT7) „ Q 8423 (7 18) 



(The other root is negative). This is again slightly larger 
than the formal exact Nf_ cr = 3/4. 

For the only physical value where there is an IR fixed 
point, Nf = 1, we thus have 



1ir,si = g = 0.6666... 
5 

7111,21 = r = 0.5555... 
y 



lm,u = = 0.37317... for N f = 1 . (7.19) 



Again, we find the same trend as for ^ in F + F, 
namely that the value of the anomalous dimension 7 m 
evaluated at the IR fixed point decreases somewhat when 
one goes from two-loop order (or the scheme- independent 
result) to three-loop order. 



VIII. CHIRAL SUPERFIELDS IN THE 
SYMMETRIC OR ANTISYMMETRIC RANK-2 
TENSOR REPRESENTATION 

In this section we analyze the UV to IR evoluation of 
the supersymmetric SU(iV c ) theory with $i, in the 
R and R representation, where R is a symmetric or anti- 
symmetric rank-2 tensor representation, denoted S2, A2, 
respectively. Since many formulas are closely related to 
each other, it is convenient to treat these two cases to- 
gether, as the Ti representation. In each of the combined 
formulas involving a ± or =p sign, the upper and lower 
signs apply to the S2 and A- 2 representations, respec- 
tively. For N c = 2, the S2 representation is the adjoint 
representation, which has already been discussed. Thus, 
for the S2 representation, the distinct cases begin with 
N c > 3. For the A2 case, -/V c is implicitly taken to be 
N c > 3, since this representation is the singlet if N c = 2. 
Further, note that the A2 representation with iV c = 3 is 
equivalent to the conjugate fundamental representation 
so, with our vectorlike content of chiral superfields, this 
reduces to the case 3>i , $i in the F + F representation al- 
ready covered above. Thus, the A2 cases that are distinct 
have iV c > 4. We focus here on jir^si and jiR2t' 



A. ft Function 



We have 



61 = 37V C - {N c ± 2)N f 



(8.1) 
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and 



37V c 2 (i_ A r /)T 8(7V c -7V- 1 )7V / 



(8.2) 



The expression for 63 is similarly obtained in a straight- 
forward manner from the general result (|2.5[) . The re- 
quirement of asymptotic freedom requires b\ > 0, i.e., 



N f < N f . blz = N f 



3N C 
N r ±2 



(8.3) 



For the S 2 representation, Nf max increases monotoni- 
cally from 3/2 for N c — 2 to 3 as N c — > 00, while for the 
A 2 representation, Nf^ max decreases monotonically from 
9 for N c = 3 to 3 as N c -> 00. 

The exact result for the lower boundary of the IR con- 
formal phase is 



Nf,, 



N 



f,max 



3N n 



2{N C ± 2) 



(8.4) 



For the S2 representation, Nf tCr increases monotonically 
from 3/4 for N c = 2 to 3/2 as N c -> 00, while for the A 2 
representation, Nf cr decreases monotonically from 9/2 
for N c = 3 to 3/2 as N c — > 00. Thus, the IR conformal 
phase exists for 



3A n 



2(N C ± 2) 
The coefficient 62 = for 

Nf = Nf lb2z 



<N f < 



3A C 
JV B ±2 



3A 2 



3A 2 ± 8(N C - Nc ) 



.5) 



(8.6) 



This is always smaller than Nf iCr for the S 2 representa- 
tion, so that &2 has fixed (negative) sign in the IR con- 
formal phase in this case. For the A 2 representation, if 



Nf < 1 + = 3.236.., then N ftb2z > N flCr , while if 
3.22 < Nf < Nf^ max , then Nf tb2z < Nf^ cr . Hence, the 
only physical case where N b2z > Nj iCr is for the inte- 
ger value N c = 3, where the A 2 + A 2 representation is 
equivalent to the F + F representation. 

At the two-loop level, the IR zero of f3 occurs at 
am, 2 i = -6i/62, i.e., 



OiiR,; 



2it[ZN c - (N c ± 2)N f ] 
3A 2 (A/ - 1) ± 8(A C - Nc X )N f 



(8.7) 



For this T 2 case, 



ci = 



B. Tn 



4(i V e ±2)(jV e Tl) 
N c 



(8.8) 



and 



C2 



4(A C ± 2){N C T 1)[A 2 T 2iV c - 4 + N C (N C ± 2)N f ] 
N 2 C 

The expression for c 3 is similarly obtained from the gen- 
eral result (US). 

Hence, 



7ir,SI 



2(N C ± 2)(A C T l)[3iV e - (Nc ± 2)N f ] 
3N3(N f - 1) ± 8(A 2 - l)iV c 



and 



(8.10) 



(JV C ± 2)(A C T 1)[3JV C - (JVc ± 2)A/][A 2 (3 + TV/) ± 2A C (3 + iV/) - 12][N c (N f - 1) ± 2Nj\ 
1IRM ~ [3A C 3 (A/ - 1) ± 8(A 2 - l)N f } 2 



ill) 



r 



One perturbative estimate of Nf :Cr is obtained by set- 
ting jir^si = 1 and solving for A/y- This gives 



N 



f,cr,SI 



3A C (3A 2 ± 2N C - 4) 
(A C ±2)(5A C 2 ±4A C -* 



3.12) 



Comparing these with the respective exact expressions 
for Nf tCr for S 2 and A 2 , we find 



N 



f,cr,SI 



Nf, cr = 



3N? 



(A C ±2)(5A 2 ±4A C 



(8.13) 



This difference is positive for all N c for both the S 2 and 
A 2 cases. Thus, as with the fundamental and adjoint 



representations, for these rank-2 tensor representations, 
this perturbative approach overestimates Nf tCr and hence 
underestimates the size of the IR conformal phase. A sec- 
ond perturbative estimate of Nf jCr is obtained by setting 
liR,2t = 1 and solving for Nf . The condition jiR t 2t = 1 
is a cubic equation in Nf, from which we extract the 
physically relevant root. This second method yields es- 
timates of Nf. cr that are qualitatively similar to those 
obtained with the first method with Jzr,SI- This quali- 
tative agreement between these two perturbative meth- 
ods for these rank-2 tensor representations is the same 
as what we found for the fundamental and adjoint repre- 
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sentations. 



IX. DISCUSSION AND COMPARISON WITH 
NON-SUPERSYMMETRIC SU(iV c ) GAUGE 
THEORY 

From our calculations on an S\J(N C ) gauge theory with 
M = 1 supersymmetry in this paper, we have found sev- 
eral general results. Our most detailed analyses here were 
for the cases R = F and R = Adj, with briefer studies 
of the S2 and A2 cases. It is useful to compare our re- 
sults with what we found in [l8| (see also whose 
results were in agreement with those in jl8l |) for a non- 
supersymmetric SU(iV c ) gauge theory with Nf copies of 
massless fermions in various representations. We believe 
that our findings for the supersymmetric gauge theory, 
besides being of interest in their own right, provide fur- 
ther insight into the results that we obtained previously 
for the non-supersymmetric theory. 

First, for the IR zero of f3, we find that am.u < otin,2i- 
This is the same type of shift that we showed earlier for 
the non-supersymmetric theory with the same R and sug- 
gests that the lowest-order (two-loop) perturbative calcu- 
lation of the IR fixed point gives a larger value than the 
true value. Second, we find that when one goes from the 
two-loop anomalous dimension evaluated at the two-loop 
IR zero of /3, jiR^e or the scheme-independent "fiR t si, to 
the three-loop result jiR^e, the value decreases. Again, 
this is the same trend that we found for the corresponding 
non-supersymmetric theory in [18] . Thus, as with the IR 
zero, this suggests that for both the non-supersymmetric 
and the supersymmetric theory with corresponding mat- 
ter field representation content, the lowest-order calcu- 
lation of the value of 7 m at the IR fixed point gives a 
larger value than the true value. The exact value of Nf tCr 
is not known for the non-supersymmetric theory, and an 
intensive research program has been underway for sev- 
eral years, especially using lattice measurements, to de- 
termine Nf tCr for a given iV c and R. Here we have taken 
advantage of the fact that Nf tCr is known exactly (at 
least with the level or rigor that is usual in physics) for 
the supersymmetric SU(_/V C ) theory. We have used one 
method for obtaining a perturbative estimate of Nf >cr 
here, namely to set 7 m = 1 and solve for the value of 
Nf where this occurs. With this method, we have found 
that, for a given N c and matter superfield content, the 
perturbative calculation yields a slight overestimate of 
Nf^ cr as compared with the exactly known value. This 
result a gree s with and complements the different analysis 
in Ref. [16| , which was based on an analysis of an approx- 
imate solution to the Dyson-Schwinger equation for the 
fermion matter field propagator. Our calculation of 7 m 
at the IR zero of the beta function provides some insight 
into this. Since, at least at the scheme-independent and 



two-loop level, 7 m increases throughout the IR conformal 
phase as Nf decreases from Nf tTnax , and since the lowest- 
order calculations yield a larger value of am and j m than 
the true value, it would follow that setting j m — 1 to de- 
termine the lower boundary of the IR conformal phase 
would yield a value of iV/. cr that is larger than the true 
value. One must, however, add the caveat that at the 
four-loop level for the non-supersymmetric theory and 
at the three-loop level for the corresponding supersym- 
metric theory, we have found that j m does not increase 
monotonically as Nf decreases from Nf^ max , which com- 
plicates the intepretation of the results. The deviation 
of 7j.r i 4£ from jiR t si in the non-supersymmetric theory 
was relatively small throughout much of the Nf interval 
of interest, but here we have found that the deviation 
of 7/a,3£ from ~/m,2£ (or Jm,si) is significant in the N f 
region of interest, which limits what one can infer from 
calculations of 7/^,3^. 



X. CONCLUSIONS 

In this paper we have considered an asymptotically 
free vectorial SU(-/V C ) gauge theory with J\f = 1 super- 
symmetry and Nf pairs of chiral superfields 
i = l,...,Nf, transforming according to the represen- 
tations R and R, respectively, where R includes funda- 
mental, adjoint, and symmetric and antisymmetric rank- 
2tensor representations. We have studied the evolution 
of this theory from the ultraviolet to the infrared, taking 
account of higher-loop corrections to the /3 function and 
the anomalous dimension j m . We have compared the re- 
sults obtained from the two- and three-loop calculations 
of the beta function and anomalous dimension (in the 
DR scheme) with exact results. In particular, we have 
calculated perturbative estimates for the lower boundary 
of the conformal phase, Nf >cr , by setting the scheme- 
independent and two-loop perturbative expressions for 
7 m equal to the rigorous upper bound (|4.4p . taken to be 
saturated at Nf jCr . We find that this perturbative calcu- 
lation somewhat overestimates Nf tCr as compared with 
the exact results, and thus underestimates the size of the 
IR conformal phase. Keeping in mind the caution that 
perturbative calculations become less reliable as the in- 
frared fixed point am gets larger, our results provide a 
measure of how closely perturbative calculations repro- 
duce exact results for these theories. 
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